Green's operator for Hamiltonians with Coulomb plus polynomial potentials 
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The Hamiltonian of a Coulomb plus polynomial potential on the Coulomb-Sturmian basis has an infinite sym- 
metric band-matrix structure. A band matrix can always be considered as a block-tridiagonal matrix. So, the 
corresponding Green's operator can be given as a matrix-valued continued fraction. As examples, we calculate 
the Green's operator for the Coulomb plus linear and quadratic potential problems and determine the energy 
levels. 
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I. INTRODUCTION 

Coulomb plus polynomial potentials, 

k 

v{r) = cur 1 = a-i/r + a Q + air + a 2 r 2 + • • • , (1) 

i=-l 

are often used to model various physical phenomena. The 
Coulomb potential, v(r) — a_i/r, describes the interaction 
between charged particles. A Coulomb plus linear potential, 
v(r) — ct—i/r + a\r, also known as Cornell potential, is the 
most common potential for modeling the confining quark in- 
teraction. It is also used in atomic physics to describe the Stark 
effect which occurs when the hydrogen atom is placed in an 
electric field. The two dimensional Coulomb plus quadratic 
potential, v(r) = a_ijr + a 2 r 2 , is related to the Zeeman 
effect; the hydrogen atom in magnetic field. The quartic har- 
monic oscillator potential, v(r) = a 2 r 2 + a±r A , is used in field 
theory to model the spontaneous breaking of symmetry. It is 
evident that there is a great deal of physics that depends on the 
precise knowledge of the dynamics of the Coulomb potential 
with various polynomial interactions. 

Over the years, several approaches have been developed to 
study some special cases of this problem (for a recent review 
see Ref. y|). To the best of our knowledge, no method has 
been proposed yet that could treat this problem with arbitrary 
potential strength <n and arbitrarily high power of k. 

In this work we calculate the Green's function of the non- 
relativistic quantum Hamiltonian with a Coulomb plus poly- 
nomial potential. The Green's function of polynomial poten- 
tials on the harmonic oscillator basis has been calculated be- 
fore JH [H. The use of the Coulomb-Sturmian basis will allow 
us to incorporate the Coulomb potential. 

If we know the Green's operator, then we have complete 
knowledge of the physical system. For example, the eigenval- 
ues of the Hamiltonian coincide with the poles of the Green's 
operator. The corresponding eigenvectors can be determined 
from the relation 



\i>n){i>n 



1 

2Tri 



G(z)dz, 



(2) 



where C encircles the eigenvalue E n in a counterclockwise 
direction without incorporating other poles. 

In our previous works, Refs. Hi, the Coulomb Green's 
operator was calculated in the Coulomb-Sturmian basis. In 



this basis, the operator J = z — H has an infinite sym- 
metric tridiagonal, i.e. Jacobi, or J-matrix structure, where 
z is a complex number. We have shown that the G(z) — 
(z — H)^ 1 = J -1 Green's operator can be calculated in terms 
of continued fractions. This Coulomb Green's operator was 
used to solve the Faddeev integral equations of the three-body 
Coulomb problem |6|] . 

In the Coulomb-Sturmian basis, the Hamiltonian with the 
potential (fl~|i is represented by an infinite symmetric band ma- 
trix. An infinite band matrix can always be considered as a 
block-Jacobi matrix with m x m blocks, where m is finite. 
Thus the continued fraction becomes a matrix-valued contin- 
ued fraction. 

This paper is organized as follows. In Section II we in- 
troduce the £>-dimensional Coulomb-Sturmian basis. In Sec- 
tion III we calculate the infinite band matrix representation 
of the Hamiltonian. In Section IV we derive the matrix 
continued fraction for the Green's operator. Some applica- 
tions are considered in Section V. First, to demonstrate the 
power of the method we show an analytically known case, 
the harmonic oscillator in two and three dimensions. Then 
we consider the Coulomb plus linear confinement potential 
v(r) — a_i/r + air in three dimensions and the Coulomb 
plus quadratic confinement potential, v(r) = a_i/r + a 2 r 2 , 
in two dimensions. 



II. THE COULOMB-STURMIAN BASIS 

The kinetic energy operator in D-dimension, with D > 2, 
is given by 



H 



1 / d 2 L(L + 1) 

2 \ dr 2 r 2 



(3) 



where L = I + (D - 3)/2. The Coulomb-Sturmian (CS) 
functions are the solutions of the Sturm-Liuoville problem of 
the Coulomb Hamiltonian 10] 



Ho-^j (r\n) = ~(r\n) 



(4) 



where b > is a parameter, n is the radial quantum number, 
n = 0, 1, • ■ ■ , and A = (n + L + 1)6. In coordinate space, the 



2 



CS functions are given by 



\n) 



r(n + 1) 



r(n + 2L + 2) 
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e- br (2br) L+1 L 2 n L+1 (2br) 



(5) 



where is an associated Laguerre polynomial. The CS func- 
tions form a basis. With (r"|n) = l/r(r\n), we have the or- 
thogonality 

(n\n') = (n\n') = (n\l/r\n') = 5 nn ,, (6) 
and the completeness relations 

£|n)(n| = £|n)(n| = l. (7) 



n=0 



n=0 



III. BAND-MATRIX REPRESENTATION 



We consider the Hamiltonian 



(8) 



on the CS basis. 

By utilizing the relations 

T a r oH 

and 



a+1 



(9) 



vw{-x)x a L°(x)L°,(x)dx = 1^+^+11,5 , 
o T(n+1) 

(10) 

we can easily calculate the overlap of the CS states 



i (n + L + 1 ) for n = n' , 

(n\n') = (n'\n) = { -^vM™' + 2L+ 1) for n' = n + 1 . 

for n' > n + 1 . 

(11) 

Then starting form Eq. and using Eqs. © and (fTTT >. we can 
derive the CS matrix elements of the kinetic energy operator 



%(n + L+ 1) for n' = n, 
(n\H \n') = (n'\H \n) = <| ^ n ' {n' + 2L + I) for n' = n+l, 

for n' > n + 1 . 

From Eqs. (O, ( fTTT i and ( fT2] > it follows that the Hamiltonian H = Hq + a_i/r + ciq is tridiagonal on the CS basis. 
By further utilizing Eqs. © and ( fTOb we can also derive the following matrix elements 



r(6n 2 + 2{L + l)(6n + 2L + 3)) 



for n' = n 



{n\r\n') = (n'\r\n) = < 



46 
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-^(2n' + 2L + 1) + 2L + 1) for n' = n + 1 

kV^V - !)( n ' + 2L )( n ' + 2L + 1) for n' = n + 2 



for n' > n + 2, 



and 



(n|rV) = (n'\r 2 \n) = < 



(12) 



(13) 



[(((10n + 2i + 4)(n + 2L + 3) + 9n(n - l))(n + 2£ + 2) + n(n - l)(n - 2))] for n' = n , 
_3 j( 4n / + 2L)(n' + 2L + 2) + (ri - l)(n' - 2)] yVfW + 2L + 1) for n' = n+ 1 , 

^(2n' + 2L) v /n'(n' - l)(ra' + 2L + l)(n' + 2L) for n' = n + 2 , 

for n! = n + 3 , 
for n' > n + 3. 

Oil 
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kV^V - !)(«' - 2)(n' + 2L + l)(n' + 2L)(n' + 2L - 1) 



Similarly, one can derive (n\r k \n') CS matrix elements for 
k > 2 as well. If fc is finite, the Hamiltonian is an infinite 
symmetric band matrix, if k = 1, it is a pentadiagonal, if 
k = 2 it is a septadiagonal band matrix. 



IV. MATRIX CONTINUED FRACTION 
REPRESENTATION OF THE GREEN'S OPERATOR 



The Green's operator G is formally defined by the equation 



J{z)G{z) = G{z)J{z) = 1 , (15) 



3 



where J(z) = z — H and z is a complex number. On the CS 
basis, this takes the form 



J2{i\J(m(i'\G(z)\i") = S i: 



(16) 



Now the operator J(z) has an infinite symmetric band- 
matrix structure. An infinite symmetric band matrix can al- 
ways be considered as a block-tridiagonal or block-Jacobi- 
matrix. So, Eq. ( TTo*l > looks like 
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(17) 
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where J n . n ' and G n . n ' are m x m block matrices, and the l's 
and O's are m x m unit and null matrices, respectively. 

Just knowing the NxN upper left corner of the full Green's 
matrix is sufficient to determine physical quantities. Let us 
denote the corresponding NxN upper left corner block ma- 
trices by jW, G( N ' and lW, respectively. If we multiply 



the N x 00 part of J with the 00 x N part of G we get the 
NxN block unit matrix 1^ N >. The sum, due to the block 
tridiagonal form of J, is reduced to three block terms 



Jn,n— lGji— X,n' ~^~Jn.nG n . n ' -\-J n . n J r \G n +\. n ' 



, (18) 



where n = 0, 1, ..TV and n' = 0, 1,.JV. If n < N, only 
terms from G^ are participating in the sum. For n — N an 
extra block matrix Gn+i.u' outside the truncated subspace is 
needed: 



JN,N-l(*N-l,n' + Jn,nGn,u' + Jn,N+iGn+1,- 

We can formally eliminate this block by writing 



ON,: 



(19) 



Jn,N-iGn-1,7i' + 

[Jn.N + >/A',A'+lGjv+l,ri'(Gjv,n') _1 ] ^JV .n' $N,m • 

(20) 

This formal elimination of the elements outside of G^ 
amounts to modifying the Jn.n block of JW. 

We can calculate Gn+ 1 {Gn.u' ) ~ 1 from another relation: 

^JV+l,ivGjv,Ti'+^iV+l,iV+lGAr+l,n'+'^JV+l,iV+2Gjv+2,ra' = 0. 

(21) 

By introducing the notation 

Cn+1 = —GN+l,n'{GN,n') (Jn+1,n) 1 , (22) 

Eq. (T% can be rearranged as 

Cn+i = (Jn+i,n+i — Jn+i.n+2Cn+2Jn+2.n+i) 1 ■ 

(23) 

A repeated application of this relation results in a continued 
fraction with block matrices 



C 



N+l 



(Jn 



+1.JV+1 — Jn+1,N+2 (Jn+2,N+2 ~ </jV+2,iV+3('^Ar+3,iV'+3 

I 



'J. 



N+3,N+2j 



'N+2.N+1 



(24) 



This matrix continued fraction does not depend on the in- 
dex n' and the correction term to Jn.n is the same for all n'. 
Therefore, we can write Eq. d20b in the form 

(Jij^ —Si.NSj.NJN,N + lCN+lJN+l,N)G < - N ^ — l^ N \ (25) 

i.e. the modified NxN block-Jacobi matrix is the inverse of 

q( n ) 

(G^) 1 = J t j — Si.N^j,NJN.N + lCN+lJN+l.N ■ (26) 

The numerical evaluation of matrix continued fractions is 
very similar to those of ordinary continued fractions. In back- 
ward evaluation we start at some K > N term, neglect the 
higher terms, and evaluate ( f24b from the inside out. If a new 
approximant with larger K is needed, we have to start the 



whole process again. On the other hand, the backward evalu- 
ation is simple and provides very accurate results. 

V. EXAMPLES 

To demonstrate the power of this method, we take first the 
harmonic oscillator 

H = H + 1/2 uj 2 r 2 (27) 

in two and three dimensions. This Hamiltonian has a septa- 
diagonal structure on the CS basis, which can be considered 
as a block-Jacobi matrix with 3x3 blocks. In our numeri- 
cal example we take to = 1 and roll up the continued fraction 
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FIG. 1 : Eigenvalues of the three-dimensional (full line) and the two- 
dimensional (broken line) harmonic oscillator as a function of the CS 
basis parameter b. 



up to the first block and calculate the determinant of a 3 x 3 
matrix. Figure 1 shows the poles of G°(z) as function of the 
CS parameter b. We took I = 0, which implies L = for 
D = 3 and L = —1/2 for D = 2, respectively. It can be seen 
that even a 3 x 3 Green's matrix provides all the eigenvalues of 
(127) . They agree with the exact results E n = uj(2n + L + 3/2) 
up to machine accuracy and the results are independent of the 
choice for the parameter b of the CS basis. 

Next, we consider the Coulomb plus linear potential in 
three dimensions 

H = H a + Z/r + ar, (28) 
and the Coulomb plus quadratic potential in two dimensions 

H = Ho + Z/r + 1/2 uj 2 r 2 . (29) 

The Hamiltonian (f28b is pentadiagonal, and can be consid- 
ered as a block-Jacobi matrix with 2x2 blocks, while (l29l 
is septadiagonal, like in the harmonic oscillator case. For the 
numerical values we take Z = — 1, a = 1, and u = 1. The 
lowest 20 eigenvalues are given in Table I. We observed a sim- 
ilar stability with respect to varying b as we did in the case of 
the harmonic oscillator. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we have shown that in the £>-dimensional 
Coulomb-Sturmian basis the non-relativistic £>-dimensional 



Hamiltonian with Coulomb plus polynomial potential has a 
band-matrix structure. A band matrix can always be consid- 
ered as a block-Jacobi matrix, and, consequently, the Green's 
matrix can be constructed in terms of matrix continued frac- 
tions. A numerically converged matrix continued fraction 
gives a numerically exact Green's operator, which even on 
a very small basis provides the exact spectrum. We have 
demonstrated the power of the method in the case of harmonic 
oscillator and obtained the exact spectrum. 

TABLE I: Eigenvalues of Hamiltonians with Coulomb plus linear 
potential in three dimensions and Coulomb plus quadratic potential 
in two dimensions. The potential parameters are Z = — 1, a = 1 
and u = 1. 



n 


Eq. J28t 


Eq. (|29) 


i 
i 


V.J I lylljj 1V01 


- 1 . o / 4 5y\jj t 


2 


2.450162895052 


1.576895542024 


3 


3.756905691262 


3.828388290161 


4 


4.855671243373 


5.963137645126 


5 


5.836029886654 


8.052626115348 


6 


6.736620996511 


10.11839697526 


7 


7.578378030294 


12.16972896261 


8 


8.374205689360 


14.21142722055 


9 


9.132754730978 


16.24628453060 


10 


9.860176266906 


18.27605894134 


11 


10.56103960914 


20.30192413905 


12 


11.23885563715 


22.32469992791 


13 


11.89639544211 


24.34497987508 


14 


12.53589461658 


26.36320650647 


15 


13.15918982353 


28.37971786276 


16 


13.76781330561 


30.39477752867 


17 


14.36306021727 


32.40859467947 


18 


14.94603779901 


34.42133786062 


19 


15.51770206715 


36.43314470188 


20 


16.07888570444 


38.44412891767 



As examples, we studied the Coulomb plus linear confine- 
ment in three dimensions and the Coulomb plus quadratic con- 
finement in two dimensions. The first is related to the quark 
confinement and the Stark effect, while the second one is re- 
lated to the Zemann effect. The exact knowledge of these 
Green's operators may facilitate the use of integral equations 
to describe quantum processes in external fields. 
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